Nondissipative Spin Hall Effect via Quantized Edge Transport 
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The spin Hall effect in a two-dimensional electron system on honeycomb lattice with both intrinsic 
and Rashba spin-orbit couplings is studied numerically. Integer quantized spin Hall conductance is 
obtained at zero Rashba coupling limit when electron Fermi energy lies in the energy gap created 
by the intrinsic spin-orbit coupling, in agreement with recent theoretical prediction. While nonzero 
Rashba coupling destroys electron spin conservation, the spin Hall conductance is found to remain 
near the quantized value, being insensitive to disorder scattering, until the energy gap collapses 
with increasing the Rashba coupling. We further show that the charge transport through counter- 
propagating spin-polarized edge channels is well quantized, which is associated with a topological 
invariant of the system. 

PACS numbers: 72.10.-d, 72.25.-b, 71.70.Ej, 73.43.Cd 



The proposals of intrinsic spin Hall effect (SHE) in a 
Luttinger spin-orbit (SO) coupled three-dimensional p- 
doped semiconductor Q and in a Rashba SO coupled 
two-dimensional electron system (2DES) |2jhave stimu- 
lated many subsequent research activities @, 0, IE IE S 
IE IE EE El 111 EE 13 EE EE H 113 The SHE may po- 
tentially provide a purely electrical means to manipulate 
electron spins without use of ferromagnetic materials or 
a magnetic field. The SHE in these systems is dissipa- 
tive because of nonzero longitudinal conductance [tj and 
exhibits nonuniversal behavior in the presence of disor- 
der [EIEEEE1E1EEE1EEEEE3> whi ch is naturally 
distinct from the conventional integer quantum Hall ef- 
fect (IQHE). In particular, it is found |fj||jll6| that the 
bulk SHE in two-dimensional Rashba model may be de- 
stroyed by any weak disorder in infinite samples. It is 
of both fundamental and practical interest to search for 
nondissipative SHE with universal properties similar to 
the IQHE, in light of that IQHE can exist in nature in 
the absence of magnetic field, as first predicted by Hal- 
dane EE]- 

A class of band insulators with SO coupling were sug- 
gested as possible candidates for nondissipative SHE 
Interestingly, Kane and Mele proposed |2CJ that the in- 
trinsic SO coupling in single-layer graphene films may 
give rise to an integer quantized SHE (IQSHE). The in- 
trinsic SO coupling conserves electron spin s z . The in- 
dependent subsystems of two spin directions <r =f and 1 
are each equivalent to Haldene's spinless IQHE model [n| 
on honeycomb lattice without magnetic field. They con- 
tribute quantized Hall conductances e 2 /h and —e 2 /h, re- 
spectively, when the electron Fermi energy lies inside the 
energy gap created by the SO coupling. While the charge 
Hall conductances cancel out, the total spin Hall conduc- 
tance (SHC) is quantized to ij s h = 2 in units of e/4-7r. 
We recall that each subsystem can be classified by an 



integer Chern number EE> which equals the Hall con- 
ductance of the subsystem in units of e 2 /h, and is con- 
served without spin-mixing interactions. Upon coupling 
the two subsystems, only the total Chern number (as 
a well-known topological invariant) is conserved , which 
is trivially zero as the total Hall conductance vanishes. 
Therefore, the conservation of electron s z appears to be 
important to the IQSHE. It is unclear whether the topo- 
logical SHE could survive if electron s z conservation is 
destroyed, e.g., by the Rashba SO coupling, which usu- 
ally exists in a 2DES due to asymmetry in the confining 
potential. Furthermore, disorder effect in the class of in- 
sulating SHE systems has not been studied so far. It is 
highly desirable to investigate these important issues. 

In this Letter, the SHE in the 2D honeycomb lattice 
model including the intrinsic and Rashba SO couplings is 
studied numerically. By using the multi-probe Landauer- 
Biittiker formula, we show that the SHC remains nearly 
quantized in the presence of finite Rashba coupling and 
disorder scattering until the energy gap collapses. We 
further show that the charge transport through spin- 
polarized edge channels is well quantized even for nonzero 
Rashba coupling, which accounts for the robustness of 
the SHE. The SHC in samples with close boundary con- 
ditions is also calculated by using the Kubo formula, 
whereby the SHE is shown to be a stable bulk effect in- 
stead of a boundary effect. Our work provides the first 
numerical demonstration of robust nondissipative SHE in 
a spin nonconservative 2DES in the presence of disorder. 
The nontrivial topological origin of this nondissipative 
transport regime is also discussed. 

The Hamiltonian for a 2DES on a honeycomb lattice 
can be written as EE EE E3 
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FIG. 1: The four-probe spin Hall bar setup used for calculat- 
ing the SHC. Filled circles represent the sites in the sample, 
and opened circles stand for the sites in the leads. Ui is the 
electrical voltage in lead I. 

+ iV R ^ c\e z ■ (<x x dy)cj + ^ e l c l t c l , (1) 

(ij) ' 

where cj = (cL, c|,) are electron creation operators, and 
er are the Pauli matrices. The first term is the usual near- 
est neighbor hopping term. The second term is the intrin- 
sic SO coupling allowed by the symmetries of the honey- 
comb lattice 0, |2(1 with i and j as two next nearest 
neighbor sites, where k is the only common nearest neigh- 
bor of i and j, and is a vector pointing from k to i. 
The third term is the Rashba SO coupling with e z a unit 
vector in the z direction. The last term describes non- 
magnetic disorder, where ti is a random on-site potential 
uniformly distributed in the interval [— W/2, W/2]. The 
distance between nearest neighbor sites is taken to be 
unity. We mention that honeycomb lattice may be re- 
alized in other microstructures as well as in single-layer 
graphene films \22\. For example, in a triangular an- 
tidot lattice created at a semiconductor heterointerface 
by using artificial periodic repulsive potentials, the elec- 
trons can be restricted into the region of a honeycomb 
sublattice 23] . If we switch off the Rashba coupling by 
setting V R = 0, Eq. JQ) reduces to a two-component Hal- 
dane's model [19j . which is expected to display a S H = 2 
IQSHE HI- 

We consider a four-probe spin Hall setup as illustrated 
in Fig. 1, where a rectangular honeycomb lattice sample 
is connected with four ideal semi-infinite leads. The ac- 
tual system used in our calculations will have the same 
aspect ratios as Fig. 1 but enlarged sizes. To specify 
the system size, the sample is divided into L y horizon- 
tal chains with L x sites in each chain, as indicated by 
the dotted lines in Fig. 1. The total number of sites in 
the sample is denoted as iV = L x x L y . For simplicity, 
the leads are assumed to have a square lattice structure 
with only nearest neighbor hopping t. The spin currents 
are well defined in the leads, where no SO interactions 
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FIG. 2: (a) The curves with filled areas are electron DOS 
(bottom axis) in a clean bulk sample versus electron energy 
(left axis) for four different Rashba coupling strengths Vr. 
The dashed lines represent the two edges of the energy gap 
as functions of the Rashba coupling strength V R (top axis), 
(b) Four-probe SHC versus electron Fermi energy for some V R 
calculated on a N = 129 x 64 sample. For both (a) and (b), 
W = and V so = O.li. 

exist. The SHC is given by twice the ratio of the spin 
current in lead 3 to the voltage drop V between leads 
and 1. Here, a factor 2 is used to properly eliminate the 
effect of the contact resistances between the leads and 
the edge channels in the four-probe setup [2(J. The lin- 
ear SHC is calculated exactly by using the multi-probe 
Landaucr-Bfittiker formula EE 

SH = £T,-( T so'-T 3 7') , (2) 

aa' 

where TfiF is the spin-resolved electron transmission co- 
efficient from spin a' channels in lead V to spin a channels 
in lead I at the Fermi energy E. 

In Fig. 2a, we show the electron density of states (DOS) 
calculated in the momentum space for a large clean bulk 
sample. The SHE within the energy gap is of our main 
interest. For a sample at half filling, such as an undoped 
graphene film, the presence of weak disorder can pin the 
electron Fermi energy inside the gap. We note that the 
actual electron DOS of the sample in the four-probe setup 
may be slightly modified from that shown in Fig. 2a be- 
cause of the connection with the leads. In Fig. 2b, the 
calculated SHC is shown as a function of electron Fermi 
energy E for sample size N = 129 x 64 and several dif- 
ferent strengths of the Rashba coupling V R . At V R = 0, 
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FIG. 3: Four-probe SHC a S H and total transmission coeffi- 
cient T30 for Vso = O.lt, E = OAt and three different sample 
sizes, (a) to (c): o s h versus disorder strength W for V R = 0, 
O.lt and 0.2t, respectively, (d) to (f): corresponding total 
transmission coefficient T30 from lead to lead 3 versus W . 
Here, disorder average is performed over 1000 random config- 
urations. 



the energy gap in the DOS for a clean bulk sample is 
between — 0.52* to 0.52*. As expected, the SHC in Fig. 
2b is well quantized to integer 2 in the main region of the 
gap. As V R increases to O.lt, the gap shrinks to — 0.23i to 

0. 511. The SHC within the gap deviates from the quan- 
tized value slightly, showing very small fluctuations with 
E, in contrast to the strong fluctuations outside the gap. 
The same feature is observed for V R = 0.2t and 0.3*. 

The effect of disorder is studied by fixing electron 
Fermi energy at E = OAt. For V R < 0.3* considered 
below, E — OAt is always inside the energy gap of a pure 
bulk system. The SHC calculated for three different V R 
are plotted in Fig. 3a-3c, respectively, as functions of W. 
The results for three different sample sizes (N = 65 x 32, 
97 x 48 and 129 x 64) are shown together for comparison. 
In Fig. 3d-3f, the corresponding total transmission coeffi- 
cient T 30 = J2aa> T io' is plotted. All T w for neighboring 
leads I and V are equal after disorder average by symme- 
try. They characterize the charge transport between the 
leads. 

Let us first consider the case of zero Rashba cou- 
pling V R = 0, where the SHC is integer quantized at 
zero disorder. According to Fig. 3a, the IQSHE per- 
sists for a range of disorder strength < W < 1.2*. 
In the same range, T30 in Fig. 3d is also quantized to 

1. This is not surprising because at V R = the sub- 
systems of the two spin directions are two decoupled 
IQHE systems. Our result is consistent with electron 
transport through fully spin-polarized edge channels with 
spin-dependent chirality. In the IQSHE regime, our cal- 



culation yields T. 
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1, represent- 
ing a left-moving edge mode in a =f subsystem, and 
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1, corresponding to a right- 
moving edge mode in a =\ subsystem. All the other spin- 
resolved transmission coefficients vanish. Strong disorder 
W > 1.2* destroys the quantizations of the SHC and 
transmission coefficients. On the strong disorder side, 
T30 increases rather than decreases with increasing W, 
which signals the collapse of the bulk mobility gap. It is 
verified but not shown here that, with further increasing 
W, all the transmission coefficients eventually decrease 
to zero because of electron localization. 

Next we look at how the SHC and charge transport 
evolve with disorder at nonzero Rashba coupling. Re- 
markably, we see from Fig. 3e and 3f that T30 is still 
well quantized within a relatively small range of W, 
indicating that the edge modes remain robust. T^o(= 
T13 = T21 = T02) = 1 relates to the left-moving mode, 
and T 03 (= T 31 = T 12 = T 20 ) = 1 relates to the right- 
moving mode. However, the spin-resolved transmission 
coefficients are no longer quantized. For example, at 
V R = 0.2t, W = 0.5* and N = 129 x 64, we have 



T 3 T T = 0.960, = 0.028, 
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suggesting that the edge modes become partially spin- 
polarized. Nonetheless, as long as the charge transport 
is quantized, the SHC stays near the quantized value 
(o s h — 1-95 for the above parameters), and is robust as 
it is insensitive to disorder W and independent of sample 
size N, as seen from Fig. 3b and 3c. With further in- 
creasing W, T30 deviates from the quantized value, and 
the SHC also decreases rapidly, the system undergoing a 
transition to a dissipativc transport regime. 

We have observed that in the presence of not too strong 
Rashba coupling and disorder, the SHE remains robust 
and nearly integer quantized. Similarly to the IQHE, 
while the effective current-carrying states are edge states 
in open-boundary systems, the nearly quantized SHE is 
essentially a stable bulk effect insensitive to boundary 
conditions or local Hamiltonian at the boundary. It is of 
interest to demonstrate this point directly, especially, in 
view of the sensitivity to boundary conditions of the SHE 
in other metallic models [12j, ll-l |l4j, |l7j • We consider a 
square sample without leads. Periodic boundary condi- 
tions are imposed in both the x and y directions. This 
close-boundary system has translational invariance in the 
absence of disorder. The Kubo formula is conveniently 
used to calculate the SHC Q 'hH by exact diagonalization 
of the system Hamiltonian |lq|. Standard spin current 
operator 0] J* s = {s z v y + v y s z )/2 is adopted, where v y 
is the electron velocity operator in the y direction. 

In Fig. 4, the dashed line is the calculated <j s h for the 
ideal case of zero disorder W = and Rashba coupling 
V R = for system size N — 64 x 64. <j s h is well quan- 
tized to 2 in units of e/47r within the energy gap (—0.52* 
to 0.52*). The lines with symbols are u s h at nonzero 
W = t and V R = O.lt for four different sample sizes from 
N = 24 x 24 to 64 x 64. In comparison with the ideal case, 
the SHC in the gap (—0.23* to 0.51* for a pure sample) 



4 



2.5 



2.0 - 




V R = 0.1t, w = t 




N = 24 x 24 % K 
N = 36x36 * x - 1 
N = 48 x 48 
N = 64 x 64 




-0.5 0.0 0.5 

Electron Fermi energy E (t) 



FIG. 4: SHC calculated from the Kubo formula as a function 
of electron Fermi energy E for Vso ~ O.lt. The dashed line is 
a S H in a clean sample (W = 0) with V R = and N = 64 x 64. 
The lines with symbols are <t s h at W = t and V R = Q.lt for 
four different sample sizes. Here, o s h is averaged over 200 
disorder configurations for TV = 64 x 64 and 1000 disorder 
configurations for the smaller samples. 



is not well quantized because of nonzero V R . However, 
the SHC is very close to the quantized value even in the 
presence of disorder W = t. The four lines collapse in 
the middle region, an indication that the nearly quan- 
tized SHC does not change with increasing sample size, 
and is thus expected to persist in the thermodynamic 
limit. We have seen that the SHC in this close-boundary 
system behaves similarly to that in the four-probe setup, 
provided that the electron Fermi energy lies inside the 
energy gap. A nonessential discrepancy is observed that 
the nearly quantized SHC at V R > for the former sys- 
tem is a little greater than 2 (<j s h — 2.07 in Fig. 4) and 
that for the latter system is smaller than 2. We believe 
that this discrepancy is caused by the definition of bulk 
spin current Jy S , which is not conservative [g and hence 
is not completely equivalent to the spin current measured 
in leads. 

The robust SHE and the quantized charge transfer 
through edge channels are associated with the nontriv- 
ial topological properties of the honeycomb lattice model 
with SO couplings. In the absence of the Rashba cou- 
pling, we have two decoupled subsystems of spin a =f 
and I, and each exhibits an integer quantized Hall con- 
ductance, e 2 /h and — e 2 /h. Each subsystem can be classi- 
fied by an integer Chern number, with C-\ — —C^ = 1. In 
an open system, there will be two decoupled chiral edge 
modes moving in opposite directions along the bound- 
ary. This picture is substantially altered by the Rashba 
coupling when the mirror-plane symmetry is destroyed. 
There are no longer two "types" of electrons due to spin- 
mixing effect of the Rashba term. However, while the to- 
tal Chern number vanishes, the opposite nonzero Chern 
numbers in the coupled system can not annihilate each 
other, as a consequence of "parity anomaly" in the de- 
coupled limit ■ They coexist and lead to a new topo- 



logical invariant, which manifests as a pair of edge modes 
with partial spin polarizations, as indicated by the nu- 
merical results. These edge modes are robust in the pres- 
ence of disorder, until the energy gap collapses, where 
the low-energy states merge with their high-energy parity 
partners [l^. It is worth stressing that the SHC itself is 
not a topological invariant, which decreases continuously 
with increasing the strength of the Rashba coupling as 
the edge states become less spin-polarized. Mathemat- 
ical description of the new topological invariant will be 
reported elsewhere. 

Note added: As we are finishing this paper, it is inter- 
esting to notice that in a couple of recent preprints [25) . 
different models for IQSHE are proposed and studied for 
pure systems. Our paper represents the first numerical 
work on the characterization of the SHE in this class of 
models in the presence of random disorder and coupling 
between different topological subsystems. 
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